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O ' Abstract 

On . 

Q\ . A special case of metric-affine gauge theory of gravity (MAG) is equivalent to 

q \ general relativity with Proca matter as source. We study in detail a corresponding 

numeric solution of the Reissner-Nordstrom type. It is static, spherically symmetric, 
and of electric type. In particular, this solution has no horizon, so it has a naked 

rr| singularity as its origin. 
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1 Introduction 



We will present a numeric solution of the Einstein-Proca theory with motivation that 
this theory is equivalent to a special case of the metric-afhne gauge theory of gravity 
(MAG) [4]. This special case is the so-called triplet ansatz [6] in which, roughly, the 
lagrangian includes square terms of nonmetricity but also square terms of the derivative of 
nonmetricity (stemming from curvature squares). It is a special feature of this ansatz that 
the nonmetricity (as well as torsion) may be expressed by a 1-form, which, because of this 
type of lagrangian, represents a massive, propagating 1-form field, i.e. perfectly equivalent 
to a Proca field. Obukhov et al. [6] gave an exact proof of this equivalence. 

The outline of this paper is rather straightforward: First, we briefly present the triplet 
ansatz and Obukhov's theorem. Then we review the Einstein-Proca theory, set up the 
lagrangian, derive the field equations, and simplify them for a spherically symmetric ansatz. 
After addressing the problem of integration constants and dimensions we can perform the 
numerical integration. We also present a power series expansion of the solution around the 
origin. The most important features of our solution are summarized at the end. 



MAG and the triplet ansatz 

^ o 

Given the curvature R^ and its 6+5 irreducible pieces ^R^ and ^R a ^, the torsion T a 
and its 3 irreducible pieces ^T a , and the nonmetricity Q a /3 and its 4 irreducible pieces 
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^Qa/3 ( see [4] appendix B), we may write the general MAG lagrangian as 
1 



£mAG ~ 7T~ 

2k 



— a R a P A r] a p — 2\r] a la Hilbert-Einstein 

+ T a A *(ai=i,,3^'T a ) quadratic torsion 

+ Qa/3 A *{bi=i.A Q al3 ) quadratic nonmetricity 

+ h ( (3) Q a7 A #*) A "(^Q^ A tip) quad. nonm. mixed with tf c 

+ 2(c/=2..4 Q a p) A "d a A *T^ cross terms nonm. /torsion 

— — R al3 A * (ttfj=i..6 Rap) + (^j=i..5 R a p) quadratic curvature 

+ w 7 $ a A (e 7 J ^ify) + z 6 ^ A (e a J ( 2 ^) 

+ z I=7 .. 9 K A (e 7 j C- 4 ^)] . } mixed with ^ (1) 

Here, we sum over the repeated index I . This lagrangian and the presently known solutions 
have been reviewed in [3]. We have the weak and strong coupling constants 1/k and 1/g, 
the cosmological constant A, and the 28 parameters a/ =0 ..3, &r=i..5, c/ =2 .. 4 , i«j=i..7, and 
2/=i..9- Note that the weak coupling constant has length dimension [1/ft] = l/^ 2 because 
it multiplies to a torsion square, torsion being the field strength of translation generators 
with dimension l/£. The strong coupling constant, though, has no length dimension. The 
work presented here is only concerned with the so-called triplet ansatz, i.e. the special case 

0=1, Wi=i.. 7 = , Zj=i.. 3 = Zr=5..9 = , z 4 ^ . (2) 

This means that we consider a general weak lagrangian (with weak coupling constant) but 
only a very restricted strong lagrangian (of curvature squares) allowing only for a square 

of the dilation curvature ^R a p '.— -g a pR~p. Qualitatively, the lagrangian with these 

constraints may be displayed as 

C ~ A + R + (T + Q) 2 + (MR) 2 . (3) 

Here, R, T, and Q denote just some terms linear in the curvature, torsion, and nonmetricity, 
respectively. On this qualitative level, the result of Obukhov et al. [6] is the following: 
Effectively, the curvature R may be considered Riemannian, T and Q may be replaced by 

O 

a 1-form 0, i.e. T + Q ~ 0, and ^R is similar to d(f>. Hence, (3) reads generically 

£ ~ A + R licm + 2 + (#) 2 , (4) 

which describes an Einstein spacetime with a massive 1-form field 0, i.e. a Proca field. 

We will review the results of [6] in detail. First, one considers a special case of the MAG 
lagrangian (1) with constraint (2) by specifying the remaining parameters A, aj=o..3, 6/=i..5, 
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ci=2.a, and z±. This choice is done in [6] eq (4.1) and turns out to effectively produce a 
purely Riemannian Hilbert-Einstein lagrangian, cf. [6] (4.6). This allows to introduce a 
new variable, the effective Riemannian curvature. Then, having investigated this special 
lagrangian, they generalize it again by adding a general lagrangian (restricted by (2)) to it, 
[6] (5.1-5.5). With the aid of the effective Riemannian curvature, the field equation FIRST 
(the variation with respect to the anholonomic coframe, see [4]) reads like an Einstein 
equation with an energy-momentum source that depends on torsion and nonmetricity, [6] 
(5.11). The field equation SECOND (the variation with respect to the linear connection) 
becomes a system of differential equations for torsion and nonmetricity alone. In the 
vacuum case, where the energy-momentum S Q and the hypermomentum A a @ of matter 
vanish, these differential equations (i.e. SECOND) reduce to 

{1) T Q = , {2) T Q = ^ $ a A , (3) T a = , cf. [6] (6.2,6.6+2.3,5.20) 

{1) Q a p = , (2) Q Q/3 = , cf. [6] (5.27,6.3) 

(3) Qa/3 = ^1 (lV/3)J<£ - \9a^) , {4) Qaf3 = Ma/3 , Cf. [6] (6.5+2.7,2.8) 

cfd0 + m 2 = O. cf. [6] (6.7) (5) 

Here, the new 1-form determines the torsion and nonmetricity completely and needs to 
solve the Proca equation. The four constants m, ko, hi, and A; 2 uniquely depend on the 
parameters of the MAG lagrangian via [6] (6.8, 6.4, 5.3-5.5). We summarize 

k = 4(a 2 - 2a )(6 3 + a /8) - 3(c 3 + a ) 2 , 

h = 9/2(a 2 - 2a )(6 5 - ao) - 9(c 3 + a )(c 4 + a ) , 

k 2 = 12(6 3 + ao/8)(c 4 + a ) - 9/2(65 - a )(c 3 + o ) , 



m = -46 4 + -a + — (b 5 -a ) + — (c 4 + o ) . (6) 

Obviously, these parameters depend only on ao, a 2 , 63, 64, 65, C3, C4, and z 4 . With (5) 
and (6) we can express the energy-momentum source of torsion and nonmetricity in the 
effective Einstein equation in terms of [6] (7.3, 7.5). This energy-momentum is exactly 
the energy-momentum of the Proca 1-form 0. Hence, we finally found that the MAG 
lagrangian (1), restricted by (2), together with its field equations, is effectively equivalent 
to an Einstein-Proca lagrangian as suggested in (4). The parameter m given in (6) has 
the meaning of the mass parameter of the Proca 1-form 0. If m vanishes, the initial MAG 
theory is equivalent to the Einstein-Maxwell theory. In general, m = is equivalent to 



= 32 64 a 2 63 — 12 a 02 63 — 64 64 ao 63 — 24 63 C4 — 48 63 C4 ao — 32 64 ao — 24 64 c 3 
9 
9 



9 a 2 65 a — 6 a 2 a 2 + 9 a c 3 2 — 48 64 c 3 a + 4 64 a 2 a — 24 a 2 c 4 + 9 a 6 5 2 



-a 2 6 5 -3a c 4 +18030465-18030400 + 18030065 + 18000465. (7) 
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This equation generalizes [7] (4.2) for 65 7^ 0. Thus, the exact solution found in [7] with 
m = corresponds to an exact solution of an Einstein-Maxwell system. Here, we want to 
present a solution for m^O. 



2 The Einstein-Proca theory 

Motivated by the previous section we now concentrate on the Proca lagrangian Cp of a 
massive 1-form 0: 

C P = --d0A*d0 + -TO 2 0A*0. (8) 

First, we shortly discuss the dimension of to. We know that the Hodge-dual of a p-form is 
an (n— p)-form. Hence, whenever the Hodge-dual applies on a p-form, it has the dimension 

[*] = l n ~ 2 P. It follows that [d<p A *d<p] = [0] 2 , and A *<p = £ 2 [(j)} 2 . To be able to 
consistently add these terms in the lagrangian the dimension of the mass parameter needs 
to be [to] = l/£. 

It is straightforward to calculate the Proca field equation and the canonical energy-momentum 
of this lagrangian with the Noether-Lagrange machinery presented in [4] . The variations 
yield: 

oCp , oCp oCp ,. , , 9 , . ,„, 

w =d m) + iw = - dd ^ m *• (9) 

H =: £„ = e a \C P - (e„j#)^ - (e«J«^ 

= - [- (e„J d(f)) A *d(f) - d(j) A (e a \ *d<p) + to 2 (e Q J 0) A *0 - to 2 A (e a J *0)] 

+ (e Q J d<f>) A *d0 - TO 2 (e Q J 0) A *0 
= - [(e a \ d<p) A *# - dcf) A (e a J M0) - to 2 [(e a \ 0) A > + A (e a J »]] . (10) 

Coupled with a Riemannian background, i.e. considering a lagrangian 
C = Cp + £ffiibcrt-Einstcm) we end up with the field equations 

= —d*d(fi + to 2 *0 , Proca equation , (11) 

= G a — nT, a =: X a , Einstein equation . (12) 

For completeness we also display the contracted Bianchi identities 

= rfS Q + r^ A S^ . (13) 

Also Obukhov and Vlachynsky [5] considered this system and indeed found the same 
solution we will find. Unfortunately, they did not publish their results until recently such 
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that the author did not know about their efforts for a long time. However, the integration 
and the results are presented here in much more detail. Before we concentrate on a solution 
of this system, we consider two modified versions of the problem. 

2.1 A Proca field on flat spacetime background 

If we consider the lagrangian (8) on flat, i.e. Minkowskian background we find two simple 
solutions. The first we find by an ansatz in analogy to the static electric monopole field 
Ail — q/ r dt: We suppose to be a static, spherically symmetric, and time-like 1-form 

<f>=^ldt, (14) 

r 

where t denotes the time coordinate and r the radius in spherical coordinates. With this 
ansatz, the Proca equation (11) becomes 

= -(-< + m 2 u cl ) (15) 

r 

and can be solved by the well known Yukawa potential 

Ueijr) q . , . , 
= - exp(— mr) . (lb) 

The parameter q is called Proca charge. A second solution we find by considering an ansatz 
in analogy to the static magnetic monopole field A mag = p (1 — cos#) dip: We set 

= w mag (r) (1 — cos 6) dip . (17) 

The field equations (11) becomes 

= — — r(l -cos0)(-iC ag + m 2 M mag ) (18) 

which is, of course, also solved by u mSug = p exp(-mr), where p might be called magnetic 
Proca charge. 
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2.2 A Proca vector field 

Here we want to clarify that there is a difference between an ansatz of as a f-form and 
as a vector. Consider the analogy of the lagrangian (8) for a vector-valued 0-form <p a : 

C P = -D(f) a A *D(f) a + m 2 (f) a A *(j> a , (19) 

w- D ^ + w^ +2D ' D ^ +2m ^- (20) 

= -(e a \D<t>P) A *D<j>f> + D^ A (e a \*D<f>p) + m 2 (e a \^) A *0/3 + m 2 ^ A (e a J*^ 
+ 2{e a \D^) A *D(j> p - 2m 2 (e Q j0 /3 ) A ^ 

= (e Q j£><^) A *Z% + £>^ A (e a J*Z%) - m 2 [(e a J^) A *0/3 ~ / A (e a J*^)] 
= (e a \ D<ffi) A *£>^ + D<$P A (e Q J*D^) + mV A (e a J*<^) . (21) 

Note that we introduced the covariant derivative D = d + T^ with the Levi-Civita con- 
nection T^. In flat space, say, we find that the signs in (20) and (11) do not agree. Also, 
the energy-momentum (21) is very different from (10). From the introduction it is clear 
that we are interested in a Proca 1-form and not in a Proca vector field. 



3 An solution with electric type Proca field 

Motivated by the Yukawa solution = - cxp (~ mr ) dt f the flat Proca equation and the 
analogy to the electric monopole, we make the ansatz that is a spherically symmetric, 
static, and time-like 1-form: 

= ^Idt . (22) 

r 

The general spherically symmetric ansatz for the coframe (i.e. implicitly for the metric) 
reads 

^ = fdt ) $ i = -Ldr ) $ 2 = rd6, &* = r fin d<p , f = f(r) , g = g(r) . 

(23) 

The meaning of the function g is illustrated by the relation rj = grj^t for the volume 
element r\ = *1. Writing the Proca field as = $(r) dt, with $(r) := u(r)/r, the energy- 
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momentum (10) of this ansatz turns out to be 






V -[-{V / g) 2 - mH 2 / f 2 ] 

|[ (<S>'/gf ~ ro 2 * 2 // 2 ] 
V 



tf 2 A £ 2 = 1 [-{<$>' /gf 
^ A £ 3 = | [-{V/gf 



m 



in: 



! $ 2 // 2 ] , 
! $ 2 // 2 ] , 



FAS 1 

\r a s Q ; 



for fi 7^ i/ , 
m 2 $ 2 



/ 2 

$'\4 



</ 



$'m$\ 2 /m$\ 4 



/<? 



./ 



(24) 



Recall that *(fl fl A E") =: T^ v represent the components of the ordinary 2nd rank energy- 
momentum tensor. The non-vanishing components of equations (11), (12), and (13) read 
in simplified form 



Po 





Pi 





Eoo 





En 





Eoi 





E22 






f 2 g'(u 



ru 



2*. II 



gr(f u 



2 2 \ 

g m u) 



gu + g(ru' — u) 

-«[(« - ru') 2 + m 2 r 2 u 2 g 2 /f] + 2r 2 {f 2 - g 2 ) + 4//V 3 - 4f 2 r 3 g'/g 
= -k[(u - ru') 2 - m 2 r 2 u 2 g 2 /f 2 } + 2r 2 {f - g 2 ) + 4//V 3 



ru >f + m 2 r 2 M 2^ 2//2 ] + 2 « 7 V 4 /(/<j - gf) - Qfifg - gf) 



f'g-gf 

-ngf A [(u 

+ 2/V/W - 0/0 - 2/ W + 2/V(/«/ - 2^/0 
B : = / 3 (ra' - w) [#« + 0(7-11' - u)] - g 2 m 2 r 2 u [g(fii - uf) + u(fg - gf)] 
Bi : = fg'{u - ru') + gr{fu" - g 2 m 2 u) (25) 

Here, P, E, and B refer to the Proca, Einstein, and contracted Bianchi equations and indices 
denote the respective components. These equations also include, for later purposes, the 
case of time dependent functions u, /, and g. Here, we neglect this time dependence and 
thus discard all time derivatives. We see that Bi is equivalent to Po and that En — E o 
simplifies the Einstein equations a lot. We end up with the following ordinary differential 
equations system of second order in u and first order in / and g: 

= f 2 g'{u - ru') + grifu" - g 2 m 2 u) , (26) 

= - K [(« - ru') 2 - m 2 r 2 u 2 g 2 /f 2 ] + 2r 2 (/ 2 - g 2 ) + 4//V 3 , (27) 

= Km 2 u 2 g z + 4f 4 rg' . (28) 



E 



11 



Po 
En 
Eoo 



These equations are equivalent to [5] eqs (3.5, 3.9, 3.10) found by Obukhov and Vlachynsky. 
We did most of the calculations with the aid of the computer algebra systems Reduce and 
Maple. The respective files can be found at www.thp.uni-koeln/~mt/work/1999diplom/. 
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3.1 Preparing the numerical integration 

First we briefly discuss the dimensions of the system. We know the dimensions of the radius, 
the gravitational coupling constant, and the Proca mass, [r] = £, [«] = £ 2 , [m] = l/£, 
respectively. Hence, we get rid of all dimensions by rescaling the radius variable r — > r/y/n 
and the mass parameter m — » \[Km. In practice, i.e. when investigating the equation 
system with computer algebra, we simply put k = 1, which is equivalent to the rescaling 
but saves us from introducing new variables. Also, we can eliminate the mass parameter 
m from (26-28) by the substitution r — > mr, f — > f/m, and g — > g/m. Instead, again, we 
equivalently fix m = 1 in the following without loosing generality. 

We now concentrate on the dimensionless ordinary differential equation system (26-28). 
Parameter m is fixed and obviously not an integration constant. The equation system is 
of first order in / and g, and of second order in u. It is easy to reduce it to an ordinary 
first order differential equation system by substituting u' — > v, u" — > v' and by adding a 
fourth equation v — v! to the system. Hence, a general integration of this system leads 
to four integration constants. In the case of a numerical approach, these constant are the 
initial values for /, g, u, and u' at the starting point of integration. Here, we consider only 
two starting points r < 1 and r^ ^> 1 for integrations from zero and infinity. We denote 
the set of integration constants by (fo,go,u ,u' ) in one case and by (foo,9oo, u oo, u'^) in 
the other. The discussion above suggests that we will find a 4-parameter set of solutions. 
But this is misleading. In principle it is possible to start integration with four arbitrary 
parameters at some point f £ 1R. But then, in general, one will find a solution of the 
equations (26-28) only in a neighborhood of f. More precisely, the convergence for r — > 
and r — > oo are 2 constraints and thus we expect only a 2-parameter set of suitable 
integration constants. 

First we discuss the limit r — *■ 0. We presume that g, u, and u' are finite at zero and 
that / diverges as 1/r. This can be motivated by the fact that such an ansatz solves the 
system of equations (26-28) as we will show shortly, or by the arguments that g needs to 
be finite because of the volume element, u and u' should be finite because of the energy, 
and / should behave as 1/r because of the analogy to the Reissner-Nordstrom solution. 
Later, also our numerical integrations confirm this assumption. To investigate the limit 
r^Owe therefore insert the expansion 

u = Mi + u 2 r + u 3 r 2 + .. , 

9 = 9i + 92r + g 3 r 2 + .. , 

f = fi/r + f 2 + hr + .. (29) 

into the system (26-28). Collecting the coefficients of 0-th order in r in each of the three 
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equations we find respectively 
= h 92 «i , 
= ff 9i- 2 h i9iul » 
= /i^ 2 . (30) 

These equations are solved by g-i = and /i = U\/y2. This means that, for r — > 0, we 
find the following approximate solution of the equation system (26-28): 



/. 



q 







V2r ' 

go = c, 

<, = 6 . (31) 

Hence, for an integration from zero we are left three parameters (q,b,c). Our results in 
section 3.4 will show that such solutions with arbitrary (q, b, c) diverge at some finite radius. 
Only a fine tuning of one of the parameters (we tuned c) makes the solution converging at 
infinity. Hence, regarding only converging, global solutions, we are left with a 2-parameter 
set of solutions parameterized by (q, b). 

For the limit r^ it is natural to require the metric function to tend to the Schwarzschild 
solution with the new mass parameter M that represents the total gravitating mass. Also, 
the Proca field should tend to the vacuum solution of the Proca equation, i.e. the Yukawa 
potential specified by the new Proca charge parameter Q. In detail we would expect 



Joo 


— 


V 


1- 


2M 

• 


9oo 


— 


1 


1 




^oo 


— 


Q 


exp 


{—mr 


< 


— 




mQ 


exp(- 



-mr oo) . (32) 

This gives a 2-parameter set of solutions parameterized by (M,Q). Since we sill find 
(numerically) that all of these solutions converge for r — > we conclude that this set of 
solutions contains all global solutions. 



A power series expansion 

There exists a simple scheme to determine all coefficients in the expansion (29): Inserting 
the expansion and considering the coefficients of i-th order (i > 2) in r in each of the three 
equations (26-28), one can solve for 1^+2, gi+2, and f\ + i in terms of Mj+2, 9j+2, and / J+ i 
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with j < i. This iteration is very easily implemented in Maple and we display here the 
result after considering the equations up to 4-th order: 

1 c 2 b 
u(r) =q + br - -— r 5 , (33) 

10 q z 

9(r)=c+~r A + ~r*, (34) 

2 f 5 f 

ft \ 1 4 1 c 2 1 2g 2 c 2 -3c 4 3 

f(r) =—=- + —=—r A -= - r . (35) 

y/2 r v/2 q 6^ q 3 

The problems is that an insertion of the power expansion (29) into the equation system 
exceeds the computer's memory resources very fast. Being limited in this way, we could 
not observe an appropriate convergence behavior for large r. 



3.2 Integration from infinity for various M and Q 

We perform the numerical integration with the standard Runge-Kutta method provided by 
the computer algebra system Maple. (In detail: We used the rkf 45 method with 15 digits, 
absolute (abserr) and relative (relerr) errors 10~ 13 , and unlimited number of function 
evaluations (maxfun).) Such an integration takes only about 10 seconds. You can find all 
calculations in the Maple-file given. 

Figure 1 represents a typical solution for an integration from infinity (r^ = 60, which is 
far enough from the Schwarzschild radius rg(M)). The gravitating mass M and the Proca 
charge Q are of the same order as m — 1. We see the metric functions / and g. To get a 
better impression of the behavior of / for r^Owe add a plot of h :— rf. The asymptote 
of h for r — * oo (h ~ ry/l — 2M/r ~ r — M) crosses the r-axis at r = M. The Proca 
function u = $/r exhibits a nicely localized density. Its derivative —u' is less instructive. 
Also the energy-momentum trace of the Proca field 

t : = -*(#* A £ Q ) = m 2 $ 2 // 2 = m 2 u 2 /h 2 (36) 

(see (24)) is quite localized. We also display its spatial integral 



OG 



s(r) := 4tt / r 2 g(r)t(r) dr satisfying / #* A S a = -s(0) dt . (37) 

J r ./spatial 

Next we vary M and Q. Figure 2 displays the metric functions h = rf and f 2 . The 
most interesting point of these plots is the following. For the Schwarzschild solution, f 2 
becomes negative within the horizon and vanishes at the horizon. Looking at the energy- 
momentum of the Proca field t := m 2 $ 2 / f 2 , we can already follow that in our case f 2 
may neither vanish nor be negative, as long as $ is finite. Hence, a Proca solution with 
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12 




Figure 1: Typical solution: Integration of (26-28) for m = M = Q = 1 from r^ = 60, 
with integration constants (32) performed by Maple. We see the metric functions /, g, and 
h := r f, the Proca function u and u', and the energy-momentum trace t and its integral 
s (divided by 100). 



finite $ may not have a horizon. For a better comparison, the lower two plots in figure 2 
display f 2 for our Proca system (left plot) and the square metric function of the Reissner- 
Nordstrom solution (right plot) for varying Proca charge and electric charge, respectively. 
The Reissner-Nordstrom solution lacks a horizon as long as we choose the electric charge 
larger then 2M, i.e. the over extreme case. For smaller charges the Reissner-Nordstrom 
solution has a horizon and the square metric function becomes negative. For our Proca 
system the behavior is similar for large Q. For smaller Q though, f 2 approaches zero 
but never becomes negative. The important result that our solution has no horizon is 
consistent with the analysis of Ayon-Beato et al. [1]. They proved that a static Einstein- 
Proca solution may not have a horizon by considering a spatial integral of the Proca 
equation. They called this a no-hair theorem for static black holes in the Einstein-Proca 
theory - or in the equivalent triplet ansatz of MAG. The lack of a horizon also means that 
our solutions have no continuous limit for m — > [m = means a Reissner-Nordstrom 
solution) or Q — > (Q = generates a Schwarzschild solution)! 

In figure 3 we plot the metric function g and the Proca function u. As we vary M, the 
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13 



10 



60, 1.50, 3.00, 8.00 h(r) : m=l, M= .60, Q= .05, .20, .60, 1.50, 3.00, 8.00 




harge=l. 50,1.75, 2.00, 2.25, 2.50 



Figure 2: Integration from r^. The metric functions h — rf and f 2 are displayed for 
various M and Q. For a better comparison, we plot the Reissner-Nordstrom solution in 
the right bottom. 



metric function g seems to distinguish inside and outside regions. Inside, g takes some 
constant value within [0, 1] decreasing with increasing M. Outside, g equals 1. As we 
vary Q, we find that larger values for Q smear this boundary between inside and outside. 
Looking at the Proca function u, as we vary M, we find that the Proca field becomes 
perceptibly non-vanishing exactly within the same boundary g exhibits. Very interesting 
is the curve for M = 0.6 and Q = 0.05 in the right plots. The Proca field vanishes as r — > 
and its derivative u' becomes positive. The metric function g approaches zero within the 
boundary instead of continuously approaching a finite g(0) as it does for larger Q. This 
behavior is different indeed and belongs to region II as we will explain in the following 
section. 
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3- 



g(r) : m=l, Q= .60, M= .05, .20, .60, 1.50, 3.00, 8.00 g(r) : m=l, M= .60, Q= .05, .20, .60, 1.50, 3.00, 8.00 

4- 

3- 




Figure 3: Integration from r^. The metric function g and the Proca function u are 
displayed for various M and Q. 

3.3 Comparing internal and external parameters 



After the explicit presentation of the spherically symmetric solution of the Einstein-Proca 
system, we want to examine how the external parameters M and Q are correlated to the 
internal parameters q, b, and c. Both, Q and q, are in analogy to the Proca charge - but 
with respect to different limits r — > oo and r — > 0, respectively. How are they related? 
The computational power of Maple allows to integrate the system for a quite large array 
of values of M and Q. For this array we calculated the values of the internal parameters 
q = w(0), b = ■u'(O), and c = g(0) and display them in figure 4. 

The first two of these plots display the internal Proca charge q. One can see that for any 
M, the internal q depends approximately linear on log 10 Q: 

q = a log 10 Q + j3 , where roughly 3 < a < 4.5 . (38) 

In the white regions of the left plot, the numerical integration could not reach the requested 
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Figure 4: The field configuration q = m(0), b = u'(0), and c = g(0) at zero is displayed 
for an array [—2 < log 10 M < 1, —2 < log 10 Q < 1] of different values for the external 
parameters M and Q. 
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Figure 5: The integration from zero for arbitrary constants q, b, and c: If c is too large 
(left plot) the Proca function u diverges to +oo, whereas if c is too small (right plot) u 
diverges to — oo. Note the small difference of 0.0012 between the two values of c. 

accuracy of 10~ 13 (relative and absolute error). The next two plots in figure 4 show b and c. 
The noisy peaks in the plot of b are at the very edge to regions where the integration could 
not reach the requested accuracy. However, we observe a smooth transition to positive 
values of b. This region II belongs to small values of M and Q as the diagram at the 
bottom illustrates. Also the plot of c clearly demonstrates this edge to region II but at the 
same time exhibits the smooth transition to this region for smaller Q. We cannot completely 
exclude that this behavior is an artifact of the numerical integration. It is interesting that 
this region represents the limit Q — ► which, as we discussed above, cannot be continuous. 



3.4 Solutions with fixed internal Proca charge 



As we already mentioned, the integration from zero is quite costly. If we start integration 
with the constants (31) with arbitrary q, b, and c, the solution diverges at some finite radius. 
Figure 5 displays the two possible divergences: the Proca function u either diverges to +oo, 
if c is large, or to — oo, if c is small. Only a fine tuning of c allows to find a global solution 
by an integration from zero. Of course, these integrations exhibit the same solution as an 
integration from oo. Since this procedure is very time expensive, we had to find another way 
to produce solutions with fixed internal Proca charge by using the relation (38) between q 
and Q. We fixed q by tuning Q for given M. This may be done very quickly because if 
we have found one solution with arbitrary q, relation (38) tells us of how to approximately 
choose Q for a given value of q. Thereby we need only about five steps to fix q on a given 
value up to an accuracy of 10~ 5 . Figure 6 shows how Q has to be chosen for different M in 
order to fix q = 1. Figure 7 displays the solutions for fixed q — 1 and different M. The plot 
of h nicely demonstrates the necessary relation h(0) = q/V% we found in (31). The plot of 
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Figure 6: How to choose Q for given M in order to fix q = 1. 

f 2 , again, demonstrates that our solutions do not have a horizon. Instead, f 2 approaches 
zero but never becomes negative. Finally, the plots of g and u exhibit the localization of 
our Proca particle within a finite radius. 



4 Other approaches 

Failure of the magnetic type ansatz 

In section 2.1 we found a solution <p = p exp(— mr) (1 — cos 9) dip for the flat Proca equation 
(11). However, with the general spherically symmetric ansatz (23) for the coframe, we find 
the 12-component of the Einstein equation: 



tflAXi= /m'H-l) 
gr 6 sin 9 



(39) 



which has only trivial solutions. Hence, there exists no magnetic analogue to the previous 
solution! 



Rosen's ansatz 



In this section we have a glance on the ansatz of Rosen 
(8) of a Proca 1-form together with the ansatz 



8]. He considered the lagrangian 



w (r)e~ lujt dt + Wl (r)e- lujt dr 



(40) 
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Figure 7: The metric functions h — r f, f 2 , and g and the Proca function u for fixed 
internal Proca charge q — 1. 

With (40) and the spherically symmetric coframe (23), the Proca equation (11) reads 



0=#A#A# 



i? 1 A $ 2 A tf 3 



-iuit 



fg 

fe~ iuJt 
g 2 r 



I £2 2 2\ • / 

w i\J m — w ) — iojw 

2icuwi + g 2 m 2 rw / f 2 — 2w' 
— (w'q — iujw'^r + (w' — iuwijrg' / g 



From the 023-component we read off 



iujw 



W-i 







(41) 



(42) 



f 2 m 2 — u 2 

which is equivalent to equation [8] (17). We substitute W\ and identify w = Wq. In flat 
space, the 123-component of the Proca equation (41) becomes 

= rw" + 2w' + (u 2 — m 2 )rw , (43) 
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which is in agreement with [8] (19). Comparing with (15), we find that this equation 
is the same as the Proca equation in flat space for a Proca 1-form w(r) dt with mass 
parameter \fm 2 — uj 2 . Hence, if we set to = m, as Rosen proposes in equation [8] (33) 
(in his notation C — 1 =>- to — k), then (43) is the ordinary, massless Maxwell equation. 
Thus, in flat spacetime, such a particle has no finite extension. This raises the question of 
how to choose the initial constraints at infinity for such a numerical integration of the field 
equations. Finally, Rosen assumes that the total gravitating mass (M in our notation) 
is equal to the mass parameter m (in dimensionless units). In general, one can hardly 
compare Rosen's work with ours because he concentrates on the idea of an elementary 
particle with finite and absolute boundary existing in the Einstein-Proca theory. Thus, he 
assumes an empty (exactly Schwarzschild) space outside the particle's boundary - and not 
a space that becomes Schwarzschild asymptotically, as we did. He calculates the solution 
by continuously (not smoothly) fitting the (Einstein-Proca) fields inside to the (purely 
Einstein) fields outside. 



5 Summary 

The introduction of this chapter explained the meaning of the coupled Einstein-Proca 
theory as an effective theory of MAG and thus motivated our analysis of this theory. 
Most interesting, we found the general condition (7) for the massless case, i.e. for the 
(restricted) lagrangian being equivalent to the Einstein-Maxwell theory. Then we derived 
the field equations (11, 12) and the energy-momentum (10) of the Einstein-Proca theory 
and displayed the (electric and magnetic type) Yukawa solution in flat spacetime. For an 
electric type ansatz we discussed the numerical integration and its integration constants 
and also offered the power series expansion (35) at the origin. We also proved the failure of 
the magnetic type ansatz. Here, we collect the essential features of the numeric solution: 

(1) In figure 1 we display the typical solution of the Einstein-Proca system for the case of 
the gravitating mass M and the external Proca charge Q being of the same order as the 
mass parameter m. 

(2) Figure 2 concentrates on the behavior of the metric function /, with i?° = / dt. We 
found that our solution has no horizon, which should also be clear from the energy- 
momentum trace in (24) and is consistent with [1]. Hence, our solution has a naked singu- 
larity. The lack of a horizon also prohibits a continuous limit to the Reissner-Nordstrom 
(m — » 0) or Schwarzschild (Q — > 0) solution. 

(3) Figure 3 focuses on the shape of the Proca particle. We found some boundary which 
is sharp for small external Proca charge Q. The larger the gravitating mass M, the larger 
the extension of the Proca particle. 

(4) Figure 4 exhibits the interesting linear relation (38) between the internal Proca charge 



q and the logarithm of the external Proca charge log 10 Q. 

(5) The b- and c-plots in figure 4 and the g- and u-plots in figure 3 suggest a different 
kind of behavior for small M and Q (region II). Note that this region represents the limit 
Q — > 0. Although the transition to this behavior is smooth, we cannot completely exclude 
it to be an artifact of the numerical integration. 
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